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Abstract 

A simple Mathematica code based on the differential realization of hard-core boson operators 
for finding exact solutions of the periodic- A'^ spin- 1/2 systems with or beyond nearest neighbor in- 
teractions is proposed, which can easily be used to study general spin-1/2 interaction systems. As 
an example. The code is applied to study XXX spin-1/2 chain with nearest neighbor interaction in 
a uniform transverse field. It shows that there are [A^/2] level-crossing points in the ground state, 
where A'' is the periodic number of the system and [x] stands for the integer part of x, when the 
interaction strength and magnitude of the magnetic field satisfy certain conditions. The quantum 
phase transitional behavior in the ground state of the system in the thermodynamic limit is also 
studied. 
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As is well-known, the finite periodic spin-1/2 chain with nearest neighbor interaction in a uniform 
transverse field is exactly solvable by using either Bethe ansatz or transfer matrix techniques. [^"^^ Sim- 
ilar spin chain models have been attracted a lot of attention recently due to the fact that they may 
be potentially helpful in quantum information processing[^~^l and realizable by using quantum dots, 
optical lattice, or spin interaction systems, etc^^*^"^^' Quantum phase transitions (QPTs) and entan- 
glement in these systems are of great interest because there are intimate links between the QPTs and 
entanglement. [^'^''^^^l Though numerical Bethe ansatz solution to the problem is possible and helpful in 
the large N limit, it is too complicated and difficult to be compiled into a practical algorithm for large 
but finite N cases. More importantly, there is still in need of a simple approach to exact solutions of 
spin systems beyond nearest neighbor interactions. In this Letter, we report our exact diagonalization 
algorithm for spin-1/2 systems written in Mathematica by using differential realization of the hard-core 
boson operators. The simple code can easily be used to study general one-dimensional spin-1/2 interac- 
tion systems, such as XY or XYZ spin-1/2 chains. As an example. The code is applied to study XXX 
spin-1/2 chain with nearest neighbor interaction in a uniform transverse field, which shows that there 
are a series of level-crossing points when the interaction strength and magnitude of the magnetic field 
satisfy certain conditions similar to the situation in the XX spin chain studied in [17]. The entanglement 
measure^"'^®'^^! defined in terms of von Neumann entropy of one-body reduced density matrix is used to 
measure the multi-particle entanglement and reveal the QPTs in the system. 

By using the hard-core boson mapping: = Sf + tSf b], S~ — Sf — iS^ bi, and 5° — 
Sf b\bi — i, where S'f [j-i = x, y, z) are spin operators satisfying the SU(2) commutation relations, the 
periodic condition S^^^-^ — is assumed, bi and 6- satisfy [bi,b^j] = Sij{l — 2b^jbj), [6^,6^] = [bi,bj] = 0, 

and — (fej)^ = 0, the Hamiltonian of the XXX spin-i chain with nearest neighbor interaction in a 
uniform transverse field, for example, can then be written as 

^ /I 1 1 \ ^ / 1\ 

where J > (< 0) corresponds to the anti-ferromagnetic (ferromagnetic) case, and /i is a uniform 
transverse field. Then, by using the differential realizations for the boson operators with b] — )■ Xi, 
bi —^di, (1) can be rewritten as 



where V is an operation to project a state with a;| = for g > 2 ( i = 1, 2, • • • , n) due to the hard-core 
restriction. One can easily verify that the differential realization with such restriction is consistent to the 
commutation relations of the hard-core boson operators. 

Because the total number of bosons, k = conserved, fc-'particle' wavefunction of (2) can 

be expressed in terms of fc-th order homogenous polynomials of with 

fI^\xu---,xn)= Y1 Cil-i,^n^i.---^i., (3) 

l<il<i2<---<ifc<iV 

where cl^j^,,,^^ is the expansion coefficient, and ( is used to label different eigenstate with the same k. 
Using (2) and (3), one can establish the eigen-equation 

HxxxF^ix,, ■■■,xn)= Ei^^Fl^'^\x,, ■■■,xn) (4) 

which is a second order linear partial differential equation and can easily be solved with a Mathematica 
code.P^l It should be stated that the first projection V at the end of (2) becomes an identical operation 
since there is no xf with q > 2 occurring in (3), while the final projection V should be considered in the 
code, which can simply be realized by setting Xi = W i after the matrix elements of the Hamiltonian 
being constructed. Though only an example with N = 8 and /c = 3 for the XXX spin chain with nearest 
neighbor interaction Hamiltonian Hxxx/J with /i = is shown in [20], it is obvious that the procedure 
shown in [20] can easily be extended to more general cases, such as XY or XYZ spin-i chain models 
with or beyond nearest neighbor interaction. It can be seen from [20] that we first construct the eigen- 
equation of XXX model Hamiltonian in the .T-represcntation. Then, we can obtain the energy sub-matrix 
for any k, which can be output to other code for diagonalization. Hence, the original 2^ dimensional 
energy matrix is reduced to Nl/{N — k)\k\ dimensional submatrices. Once the eigenenergy eI^^^/J and 
the corresponding eigenvector are known after diagonalization, the final wavefunction can be 

expressed as 

|fc;C)=Ff (6l,---,6]v)|0>, (5) 

where |0) is the boson vacuum and thus the SU(2) lowest weight state with 5^10) = V ?'. 

As an example of application of the code, in the following, we study quantum phase transitional 
behavior of the finite periodic XXX spin-i chain with nearest neighbor interaction in a uniform transverse 
field. One can verify that there is no quantum phase transition for the ferromagnetic case with J < 0, 
in which the ground state of the system with J < keeps unchanged in the variation of the magnitude 
of the magnetic field. Quantum phase transition occurs only in the anti-ferromagnetic cases with J > 0, 
which will be considered in the following. In order to investigate QPT behavior of the system for J > 0, 
we set J = 1 — t and h = t with < t < 1. It is clear that the ground state of the system is in the 
ferromagnetic (unentangled) phase when t = l and in the anti-ferromagnetic long-range order (entangled) 
phase when t = Q. Therefore, t serves as the control parameter of the system. In the XXX case, in 
addition to 5o = A: — N/2, the total spin of the system S is also a good quantum number. Therefore, the 
wavefunction (5) can further be written as \So = k — N/2; S, ^), where the additional quantum number ^ 
is used to label different eigenstate with the same S and 5*0. Though one can only obtain an eigenstate 
with fixed Sq from the code, one may get information about the total spin S by acting on the total 
spin lowering operator S~ = X^iXi "^i" ^o the state. For example, the state |5o = —N/2;S = N/2,^) 
must satisfy S'ISq = -N/2;S = N/2,^) = 0, while \So = I - N/2;S = N/2 - 1,0 must satisfy 
S'~|S'o = 1 — N/2; S = N/2 — 1,^ = 0, and so on, which enables us to find the corresponding quantum 
number S for each eigenstate. For A'' odd cases, the eigenstates with So = k — N/2 and S = N/2 — k 
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for fc 7^ are doubly degenerate. In such cases, the expansion coefficients C'(^) — {C'lfij-.-i^} with ^ = 1 
and ^ = 2 obtained from the code are not orthogonal with each other. In such cases, we use the Gram- 
Schmidt orthogonahzation procedure to set C'{S, = 1) = C{$, = 1) - C{$, = 1) • = 2)(7(^ = 2) and 
keep = 2) unchanged after normalization. 

It is well known that the ground state of the anti-fcrromagnetic XXX spin chain is never degenerate 
with S = ior N even and four-fold degenerate with degeneracy equal to 2(25+1) and S* = 1/2 for N odd, 
which all correspond to i = 0. We have verified that the ground state energy of the system is related to the 
following set of eigen-energies: Eg~^_lg~^-^^{t) = E!^i^{t) for = 0, 1, ■ • • , [N/2], where [a;] stands for the 
integer part of x. It should be stated that the ground state energy at t = 1 corresponds to E^^{t), while 
that at t = corresponds to E^j^^^^\t). Hence, it is clear that there are also [N/2] + 1 different ground 
states which are mutually orthogonal with the corresponding ground state energy E!^^{t), E!^^{t), ■ ■ ■, 
E^jf^^'^\t) when the control parameter t changes from 1 to similar to the situation of XX spin-| 
chain reported in [17]. Obviously, the quantum phase transitions occurring in such cases are of the first 
order. [^''J It can be verified by the code that all levels with eigenenergy E!^^^{t) for k = 1, 2, • • • , [iV/2] are 
not degenerate for N even and <t < 1, while they are all two-fold degenerate for N odd and <t < 1. 
The degeneracy of the ground state for N odd at t = is 2(25' + 1) = 4, while the ground state for N 
odd at t = 1 is a singlet with S = N/2 and = -N/2. 



Table 1. [N/2] level-crossing points for 2 < < 12. 
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.(3) 




.(5) 


.(6) 


2 


0.666666 












3 


0.600000 












4 


0.500000 


0.666666 










5 


0.566915 


0.644004 










6 


0.499123 


0.566401 


0.666666 








7 


0.511933 


0.623396 


0.655288 








8 


0.343259 


0.570166 


0.643104 


0.666666 






9 


0.462701 


0.591992 


0.642284 


0.659828 






10 


0.297378 


0.527473 


0.614872 


0.652704 


0.666666 




11 


0.420934 


0.559842 


0.621991 


0.650981 


0.662104 




12 


0.262455 


0.490059 


0.58657 


0.634069 


0.657415 


0.666666 



The first order phase transition in the system occurs due to the ground state energy level-crossing 
of -B^i„(t) with SHi(t) for i = 0, 1, • • • , [N/2] - 1 with the corresponding critical point which 
is the root of the simple linear equation (t^'^/^'"'^) = E'+^{ti^^^^^^'^) for ?: = 0. 1, 2. ■ • ■ . [A^/2] - 1. 
There are [N/2] such level-crossing points indicating that there are [iV/2] -|- 1 different ground states 
within the control parameter range < f < 1. Fig. 1 clearly shows the ground state level-crossings in 
the entire control parameter range for TV = 2,4, 5, 6, 8, and 12 cases. It is obvious that there are [N/2] 
level-crossing points dividing the ground state into [N/2] + 1 different parts, of which each is within a 
specific t range when A'' is a finite number. With N increasing, however, these specific ranges become 
smaller and smaller, and finally tend to infinitesimal, thus the ground state level becomes a continuous 
phase before crossing to E^^^ level. Therefore, there will be only one obvious critical point when N ^ oo. 
One can verify that the critical point t^^^^^^ = 2/3 is A''-independent for N even, while it will tend to 
2/3 for N odd when N ^ oo. Nevertheless, other level-crossing point tc^ values are A/'-dependent, of 
which some examples are listed in Table 1. 

Entanglement measure in the model is one of important quantities to characterize its QPT behavior, 
and is often studied by using block-block entanglement defined in terms of von Neumann entropyl^^l or 
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by using Wootters concurrenceP^l, e. g., that shown in [22]. In the following, we use the simple measure 
proposed in [17-19] with 

1 ^ 

if all N terms in the sum are non-zero, otherwise ri{'^) = 0, where ^' stands for the ground state wavcfunc- 
tion and is the reduced density matrix with i-th spin-^ fermion only. It has been shown'^''"^^'^^'^^! 

that (6) is also suitable to measure genuine A''-body entanglement in a quantum many-body system. We 
observed that is i-independcnt for the ground state in the system for N even cases, while it becomes 

i-dependent for A'' odd cases. Hence, the entanglement measure r] for A'' even cases can be simply defined 
by the reduced von Neumann entropy for any site, while it should be calculated separately for A'' odd 
cases. Table 2 shows ground state entanglement in different t ranges for N = 2, • • • , 6, respectively, in 
which the entanglement type of the ground state in each t range is indicated. For example, the state is a 
linear combination of several GHZ-like states for A'' = 4 with < t < 0.5, while it consists of two-fold de- 
generate pair which are all linear combinations of serval W-like states for N = 5 with < t < 0.566915. 
It is clear that the ground state entanglement measure gradually increases while the control parame- 
ter t decreases, which is also characterized by the quantum numbers S and S'^. In the ferromagnetic 
(unentangled) phase, S = N/2 and 5° reaches its lowest value with 5° = —N/2, while S = Sq = 
(S = -So ^ 1/2) when t < tP for A^ even (odd), in which the spin-up and -down fermions are most 
strongly correlated in comparison to that in other phases. In the most entangled long-range order phase, 
A'' even systems are most entangled with rj = 1 which is always greater than those of the nearest N odd 
systems. Furthermore, the degeneracy is doubled at the level-crossing points t = tc ■ For A^ even cases, 
the ground state is not degenerate if the control parameter t does not at those [A^/2] level-crossing points, 
while it becomes two-fold degenerate when t = t^J^ for any j due to the level- crossing. For A'' odd cases, 
the ground state is four-fold degenerate at t = and is a singlet when t > ic'^^^'"*- Besides those two 
cases, the ground state is two-fold degenerate with 5*0 = —S = —k + N/2 for fc = 0, 1, 2, • • • , [N/2] if the 
control parameter t does not at those [N/2] level-crossing points, while it becomes four- fold degenerate 
when t = t^J^ for any j due to the level-crossing. However, these degenerate states at the level-crossing 
points are still distingTiishablc from each other by the quantum number 5* and with their difference 
A(5'°) = A(S') = ±1 and by values of the entanglement measure of the degenerate states. As a conse- 
quence, for A'' even case, the ground state is not degenerate when t = 0; it becomes two-fold degenerate 
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everywhere when the control parameter t is within the half-open interval t G (0,2/3] because the level- 
crossing points are dense everywhere in this control parameter range in the N —>■ oo limit; and finally it 
becomes not degenerate again when 2/3 < t < 1. For A'^ odd case, the ground state is four-fold degenerate 
for t being within the closed interval t G [0, 2/3] in the N ^ oo limit; and it becomes not degenerate when 
2/3 < t < 1. Nevertheless, the property of the degenerate states at t = and that within < t < 2/3 are 
different for N odd case. The four-fold degenerate states at t = come from the double occurrence of 
5=1/2, while two states from So = -S = -k + N/2 and another two from So = -S = -{k +1) + N/2 
to form the corresponding four-fold degeneracy for < t < 2/3. However, it has been proved at least for 
small N cases that GHZ- and W-type states are inequivalent under the SLOCC transformations. [^'^"^^1 
Therefore, the ground state should be classified into three phases in the thermodynamic limit for A'' even 
case under the SLOCC. These three phases are one non-degenerate entangled GHZ-type phase at t = 
with T] = 1, one two-fold degenerate entangled W-type phase with t € (0, 2/3] and < t] < 1, and one 
non- degenerate fully separable phase with t G (2/3, 1] and 77 = 0. But such QPT classification is only 
meaningful under the SLOCC. For N odd case, the situation is different. There is one four-fold degen- 
erate entangled W-type phase with t G [0,2/3] and < r; < 1, and one non-degenerate fully separable 
phase with t G (2/3, 1] and 77 = 0. 



Table 2. Ground state entanglement with each quantum phase for = 2, • • • , 6. 



N 



Entanglement type in each phase 



S=-So = l 
2 Fully separable 
(?7 = 0) 2/3 < t < 1 



S = -So=0 
Bell (r? = 1) 
< t < 2/3 



S = -So = 3/2 
Fully separable 
(?7 = 0) 0.6 < t < 1 



S = -So = 1/2 
1^ = 1) is partially separable 

iVi = 0) 
14 = 2) is a W Combination 
(??2 = 0.739447) 
< i < 0.6 



S=-So = 2 
Fully separable 
(r; = 0) 2/3 < t < 1 



S=-So = l 
W (?7 = 0.811278) 
0.5 < t < 2/3 



S = -So = 
GHZ Combination 

('7 = 1) 
< t < 0.5 



S = -So = 5/2 
Fully separable 

(77 = 0) 
0.644004 < t < 1 



S = -So = 3/2 
14 = 1) is a W Combination 

= 0.610281) 
14 = 2) is a W combination 
(772 = 0.619557) 
0.566915 < t < 0.644004 



S = -So = 1/2 
14 = 1) is a W combination 

(?7i = 0.858927) 
14 = 2) is a W combination 

(r?2 = 0.858501) 
Q<t< 0.566915 



S = —So =3 
Fully separable 
(77 = 0) 2/3 < t < 1 



S = —So = 2 
W (77 = 0.650022) 
0.566401 <t<2/3 



S = -So = l 
W Combination 
(7? = 0.918296) 
0.499123 < t < 0.566401 



S=-So = 
GHZ Combination 

('7 = 1) 
0<t< 0.499123 



In summary, a Mathematica code based on the differential realization of hard-core boson operators 
for constructing energy matrix of the periodic- A'' spin- 1/2 systems with or beyond nearest neighbor 
interactions is proposed, which can easily be used to study general spin-1/2 interaction systems, such as 
XY or XYZ spin-1/2 chains. As an example. The code is applied to study the anti-ferromagnetic XXX 
spin-1/2 chain with nearest neighbor interaction in a uniform transverse field. The study shows how 



5 



the ground state of the model evolves from the ferromagnetic phase to the anti-ferromagnctic long-range 
order phase with decreasing of the control parameter t introduced. In addition, we have shown that there 
are [N/2] level-crossing points, of which the middle part will become a continuous one in the large- A'' 
limit leading to the three-phase result in the thermodynamic limit for N even case under the SLOCC, 
while there are only one entangled W-type phase and one separable phase in the large N limit for N odd 
case. Such level-crossing should be common in other spin interaction systems in a uniform transverse 
field. 
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